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We investigate non-equilibrium behavior of driven dissipative systems, using the model presented 
in [Phys. Rev. Lett. 93, 240601 (2004)]. We solve the non-Boltzmann steady state energy distri- 
bution and the temporal evolution to it, and find its high energy tail to behave exponentially. We 
demonstrate that various measures of effective temperatures generally differ. We discuss infinite hi- 
erarchies of effective temperatures defined from moments of the non-exponential energy distribution, 
and relate them to the "configurational temperature" , measured directly from instantaneous particle 
locations without any kinetic information. We calculate the "granular temperature", characterizing 
the average energy in the system, two different "fluctuation temperatures", scaling fluctuation- 
dissipation relations, and the "entropic temperature", defined from differentiating the entropy with 
respect to energy. 
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I. INTRODUCTION 

Systems of many particles interacting dissipatively are 
far from thermodynamic equilibrium, and a general the- 
oretical description of their statistical mechanics is lack- 
ing, in contrast to systems in equilibrium, for which there 
is a well established theory. In this Paper we use exact 
solutions of a simple stochastic model in order to explore 
aspects of dissipative systems. In particular, we are inter- 
ested in steady states: their non-Boltzmann energy dis- 
tribution, the way in which a system arrives at its steady 
state, and various proposed definitions of temperature. 
Previous theoretical research on energy distributions and 
effective temperatures in various driven dissipative sys- 
tems (granular materials Q p,!!, |3Jp, IE S > foams 
[lOl ITll IT^ . and glasses |l3l Il4ll5j|) has been mostly 
numerical or approximate. Therefore the resulting dis- 
tributions and temperatures may agree with an effective 
equilibrium behavior due to numerical error. 

In systems comprised of macroscopic particles, en- 
ergy is dissipated via interactions, being transferred from 
macroscopic degrees of freedom (such as motion of parti- 
cles) into microscopic degrees of freedom (heat), and can 
not be transformed back. Continuous driving is needed 
in order to maintain such a system in a dynamic state. 
One way to model this driving is by holding the sys- 
tem in contact with a bath, or large energy reservoir. 
As in 'l6^ , we concentrate on driving mechanisms where 
this bath is in equilibrium at some temperature Tg (this 
driving bath is not necessarily in thermal equilibrium). 
A non-dissipative system driven by a thermal bath would 
reach thermodynamic equilibrium with it, where the en- 
ergy distribution is given by the exponential Boltzmann 
distribution, and the system temperature is equal to the 
bath temperature. 

The temperature of an equilibrium system is mani- 
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fested in various measurements that can be performed 
on it. A non-equilibrium system does not, a-priori, have 
a unique well-defined temperature, and each such mea- 
surement inspires the definition of a corresponding ef- 
fective temperature. For example, the entropic temper- 
ature Ts is the inverse of the derivative of entropy with 
respect to energy 0, 0, 0, in analogy with the 
definition of temperature in statistical mechanics. The 
principle of energy equipartition in thermal equilibrium 
motivates defining the granular temperature Tq as the 
average energy per degree of freedom |19|: The equi- 
librium fluctuation-dissipation theorem suggests defining 
the fluctuation temperature Tp as the ratio of fluctua- 
tion to response |23, 0| ; Recent results expressing the 
equilibrium temperature from ensemble averages of par- 
ticle locations (without any kinetic information) |23, H^l 
lead to the definition of a configurational temperature 
Tc as another measure of effective temperature for non- 
equilibrium systems 0, 12^ . All these deflnitions yield 
the same value in equilibrium. 

In steady states of driven dissipative systems all effec- 
tive temperatures are generally much smaller than Tb 
(Tb should not be confused with the actual temperature 
of the environment which is typically much lower), and, 
unlike thermal equilibrium, their values depend on the 
details of the coupling with the bath. Although these 
systems are far from equilibrium, and their energy dis- 
tributions differ signiflcantly from the Boltzmann distri- 
bution, there is evidence for coincidence of different ef- 
fective temperatures. Different Tps of the same system, 
obtained from correlations and response of different vari- 
ables, have been found in numerical experiments to co- 
incide in glasses ^[ [l4L llSj , and to coincide with Tn in 
granular gases 0,11] and with Ts in sheared foam ■ 

In this Paper we study an exactly solvable dissipative 
model, in which interactions occur randomly and redis- 
tribute en erg y stochastically between the interacting par- 
ticles. In |2y] we first introduced the model, discussed its 
similarities to granular gases, calculated all moments of 
its steady-state energy distribution, and showed that Tp 



differs from Tq . Here we present these results in more de- 
tail and provide further results on the model: We define 
the model in Sec. ^and investigate its energy distribu- 
tion in Sec. IIIII We use the energy scales defined by 
high moments of the energy distribution to show that 
the high energy tail is exponential with a decay rate cor- 
responding to the bath temperature Tb, and solve the 
temporal evolution to the steady state. Section Hvl deals 
with various effective temperatures in our model. We 
show that spontaneous fluctuations of different quanti- 
ties are scaled by different fluctuation temperatures. We 
demonstrate that for systems with a single energy scale 
and with smooth energy distributions, the entropic tem- 
perature coincides with the granular temperature, and 
then investigate cases where they differ. We relate the 
infinite hierarchies of different effective temperatures de- 
fined from the energy distribution to generalizations of 
the configurational temperature. All results are calcu- 
lated exactly in the context of our model, and in Sec. El 
we discuss the generality of these results to other driven 
dissipative systems. 



II. MODEL DEFINITION 

At its most basic, a dissipative system consists of a 
set of degrees of freedom, or modes of excitation, which 
interact among themselves and with the external environ- 
ment. When two modes interact with each other, there 
is energy exchange, with some of the energy being lost 
to the environment. External driving may be thought of 
as the injection of energy into the system's modes from 
the environment. Our model is constructed as a minimal 
model including these essential features. When two par- 
ticles (we shall refer to the modes as particles for simplic- 
ity) interact they lose some of their energy and exchange 
what remains, and when a particle interacts with the en- 
vironment it is more likely to gain energy rather than to 
lose. 

Our model consists of N particles with energies {ei}, 
with a constant interaction rate between any two par- 
ticles in the system. In every interaction two particles 
from the system are chosen at random and their ener- 
gies are summed. In the case of conservative dynamics 
(analogous to elastic collisions) this total energy is repar- 
titioned randomly between the two interacting particles 
(as in (27.] ). For dissipative dynamics with inelastic col- 
lisions, only a fraction < a < 1 of the total energy is 
repartitioned between the particles and the rest is dis- 
sipated out of the system. Thus, a is analogous to a 
restitution coefficient. The system is coupled to a heat 
bath so that it may be maintained in a nontrivial steady 
state. The interactions are shown diagrammatically in 
Fig. m and described in further detail below. 

For conservative dynamics {a = 1) an isolated sys- 
tem (i.e., not in contact with the driving bath) reaches 
thermodynamic equilibrium with the exponential Boltz- 
mann distribution for each particle's energy, p{e) = 
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FIG. 1: The possible interactions in our model: dissipative 
two particle interaction (left) and conservative system-bath 
interaction (right). 



exp(— e/T), where the temperature equals the aver- 
age energy per particle T — (e) (we measure temperature 
in units of energy and set Boltzmann's constant to one). 
Dissipative dynamics (0 < a < 1) cause energy to decay, 
therefore we drive the system by attaching it to a heat 
bath, constructed as an infinitely large system of parti- 
cles obeying the conservative dynamics described above, 
kept in equilibrium at a temperature Tb- The coupling 
of the dissipative system to the bath is through conser- 
vative interactions between a particle chosen at random 
from the system and a particle chosen at random from 
the bath (the system-bath interactions are taken as con- 
servative for simplicity, however dissipative interactions 
may as well be considered, yielding qualitatively simi- 
lar results). This contact is characterized by a coupling 
strength, < / < 1, defined as the fraction of every 
particle's interactions with the bath out of all its interac- 
tions. Unlike thermodynamic equilibrium, the dissipative 
system's steady state depends on the bath through both 
Tb and / (see Q which emphasizes the importance of 
coupling details). 

The stochastic evolution of the energy of particle i dur- 
ing an infinitesimal time step dt is hence given by 



e^{t + dt) 



Value : 
e^{t) 
za [ei{t) + ej{t)] 
z {ei(t) + es] 



Probability : 
1-Tdt 
(1 - f)Tdt 
fVdt 



(1) 



where F is the interaction rate per particle per unit time, 
j S {!,..., A^} (j 7^ i) is the index of the particle with 
which particle i may interact (chosen randomly at ev- 
ery interaction), z e [0, 1] is the fraction of repartitioned 
energy given to particle i in the interaction (chosen ran- 
domly with a uniform distribution at every interaction), 
and cb is the energy of the bath particle with which par- 
ticle i may interact, which at every interaction is chosen 
randomly from the equilibrium distribution in the bath: 
PB{eB) = Tg^ exp(-e_B/TB). 

The simplicity of our model derives from the fact that 
every particle in it is described only by its energy, as 
opposed, for example, to the 2d degrees of freedom per 
particle in a d-dimensional frictionless hard sphere gas. 
By eliminating the momentum and spatial variables and 
using only the energy, we replace the vectorial collisions 
between particles by scalar interactions. Furthermore, 
since any two particles may interact, there are no spatial 
correlations. 
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III. ENERGY DISTRIBUTION 

In this section we investigate our model's single- 
particle energy distribution p{e). We calculate from Eq. 
([T|) the temporal evolution of any moment of p{e) . This is 
used to obtain the temporal evolution of p{e) , as well as 
its form in steady state. We use the energy moments to 
define two hierarchies of energy scales, whose asymptotic 
behavior is then used to characterize the high energy tail 
of p{e). In particular, we shall show that although p(e) 
differs from the Boltzmann distribution, its high energy 
tail is exponential. 



A. Average Energy 

The most direct way to characterize the system's state 
is by the average energy per degree of freedom. It is 
solved by averaging Eq. over the stochasticity in the 
dynamics and over all particles: 



{e{t + dt)) = {1-Tdt){e{t)) + {1- f)Tdta{e{t)) 
1 



+ fTdt-[{eit))+TB]. 
Hence (e(i)} satisfies the differential equation 



2 d{e{t)) 
f dt 



(2) 



(3) 



with Ai{a, /) = 1 + (1 - /)(1 - 2a) > 0. This has the 
steady state solution 



(e) = 



Tb 



ITb ^ 

Ai 2a- 1 + 2(1 -a)//' 



(4) 



In analogy with granular materials, this is denoted as the 
granular temperature Tq = (e). It is plotted vs. a and 
/ in Fig. 12 

It is interesting to note that despite the simplicity of 
our model, it captures, at least qualitativ ely, some as- 
pects of actual driven dissipative systems [23, • For 
instance, the stationary value of (e) is always smaller 
than Tb and depends not only on the dissipation through 
the restitution coefficient a, but also on the details of the 
coupling to the bath through the coupling strength /. (e) 
coincides with Tb only in the two non-dissipative limits: 
conservative interactions (a = 1) and strong coupling to 
the bath (/ = 1). 
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FIG. 2; Contours of the ratio (e) /Tb between the average 
energy per particle in the system and the bath temperature 
vs. the restitution coefficient a and the coupling strength / 
to the bath, as given by Eq. @. 



To- However, this is not the case: we find that the en- 
ergy distribution is clearly non-exponential, in qualitative 
agreement with realistic driven dissipative systems, such 
as granular gases P, 0, 0, 13 ^-^d colloidal suspensions 

[IflMllllsa. 

The first moment of p(e) is the average energy calcu- 
lated above. Higher moments are obtained by taking the 
average of the nth power of Eq. This yields the fol- 
lowing differential equation for (e") in terms of all lower 
moments and the moments of the energy distribution in 



the bath (for which (e'g) = m\T^ 



n + 1 d(e") 



= -An{e^) 



n-l 



r dt 

+f{el), (5) 



where 



A„(a,/) = n+(l-/)(l-2a"). 



(6) 



Any initial distribution will evolve with time to the 
steady state distribution given by 



B. Energy Fluctuations 

We now consider the energy distribution around the 
average energy (e). This tests whether the system is 
equivalent to an equilibrium system at an effective tem- 
perature Tg — (e). That is, whether the effect of the 
dissipation is to modify the Boltzmann distribution only 
by changing its characteristic temperature from Tb to 



(e") = { 



n-l 

E 

m— 1 



(e") [(l-/)a"(e"-'") + /(er'")] 

+f{el)}/An.{7) 



These expressions for all energy moments are exact arbi- 
trarily far from equilibrium (for general values of a and 
/) and contain information about the entire energy dis- 
tribution. In the equilibrium limits (a = 1 or / = 1) Eq. 
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^ yields the moments of the exponential Boltzmann dis- 
tribution, for which (e") = nlTg. 

The steady state energy distribution with moments 
given by Eq. iQ is shown in Fig. 13 At low ener- 
gies (of the order of several times Tq) the distribution 
is roughly exponential with a decay rate corresponding 
to the average energy Tq- At intermediate energies the 
distribution seems to exhibit an overpopulated high en- 
ergy tail, decaying slower than exp(— e/Tc). However, 
for higher energies (e > WTq = STg for the parameters 
in Fig. 13), the distribution exhibits an exponential decay 
of the form 



p{e) exp(-e/TB). 



(8) 



The overpopulation of the high energy tail is only with 
respect to scaling the energy with the average energy in 
the system and considering energies comparable to (e), 
as is customarily done in granular gases [H. I^L II. IM Ifl. 12^ . 

Interactions with the bath dominate the high energy tail, 
since a particle is much more likely to arrive at such high 
energies due to a conservative interaction with the bath 
rather than due to a dissipative interaction within the 
system (and the fraction of very high energy particles 
in the bath is larger than in the system). In the next 
section, we use the asymptotic behavior of energy scales 
defined by the moments to show that the tail is indeed 
exponential. 
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FIG. 3: Steady state energy distribution for restitution co- 
efficient a = 0.5 and coupling strength / — 0.5 (resulting in 
To = Tb/2) obtained in a numerical simulation of the model 
(solid line). Exponential distributions at temperatures Tb 
and To are given for reference (dashed lines). 



C. High Energy Tail 



non-exponential energy distribution, infinite hierarchies 
of effective temperatures may be defined from the energy 

(n) 

distribution. We define two such hierarchies, Ti, from 



(n) 



the ratios of succeeding energy moments, and Tj^j 
scaling the moments themselves: 



by 



T 



(n) 
M 



(9a) 



(9b) 



Both reduce in equilibrium to the system's temperature 
for any n. Away from equilibrium they typically differ 
and depend on n (see Fig. ^ , and for n = 1 both reduce 
to the granular temperature Tq = (e). 

The large n behavior of T^"'' and tI"-* reflects the 
distribution's high energy tail. For an exponential tail 



p(e) ~ exp(-e/re//), T^°°^ = 
stretched exponential tail, p{e) 
T[f = ff a > 1, and T^^^ 
a power-law tail T^'^ 



T 



(oo) 
M 



For a 



exp(-ce'^), 4°°) 



(oo) 



M 



OO if a < 1; For 



oo while T^"^ is undefined. 



We now show that the energy moments in our model 
are consistent with an exponential high energy tail with 
a decay constant corresponding to the bath temperature 
Tb- Since < a < 1 for n > 1, yl„ « n [see Eq. ©]. 
Furthermore, for < / < 1, since < a < 1, we have 
(1 - /)a"(e"-"} < /(e^""}, and Eq. flj) reduces to 



(e") 



-1 



/(r.-l)!}^^AT— , 
^-^ ml 

m=0 



(10) 



where the term f{e%) in Eq. lO has been incorporated 
as the m — term in the summation in Eq. (|10|l . The 
solution for large n may be approximated by taking the 
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The exponential energy distribution of a system in 
equilibrium contains a single energy scale - the tempera- 
ture of that system. For a non-equilibrium system, with a 



FIG. 4: (Color online) The hierarchies of eff'ective temper- 



atures T^"' and T^^' for several model parameters, obtained 
by substituting Eq. Q in Eq. @. 
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continuum limit, where the sum in Eq. H10() transform 
to an integral, and the resulting equation may be solved 
to yield 

(e") « CnlTJ^nf-^ (11) 

with C a dimensionless constant independent of n. 
Therefore, in the large n limit T^"'' and T^^^ both con- 
verge to Tb, 

T^^^^Tsf^y ' ^Tb (12a) 



t'^^^ ^TB{Cnf-^Y''' ^Tb (12b) 

as can also be seen in Fig. 0] This supports the observa- 
tion that the high energy tail behaves as exp(— e/Te). 

It is intriguing to speculate on the generality of this re- 
sult, that for very high energies the distribution behaves 
as an equilibrium distribution with a temperature equal 
to the bath temperature Tb (see also where similar 
results have been found for a granular gas driven by an 
ideal gas heat bath). 

D. Approach to Steady State 

We solve Eq. (O recursively with n, and find that the 
time-dependent solution is of the form 

where (e") are the steady state moments given by Eq. 
{TI), and {Cn,m] are constants depending on the initial 
distribution. Am is discrete and increases monotonically 
with m, therefore the slowest exponential decay with 
time, exp (— Airt/2), dominates the long time behavior 
of all moments. When scaling high moments to units of 
energy, one has 

Tl"^-(e"(t))^/"-exp(^-^f). (14) 

Thus the high energy components of the distribution ap- 
proach their steady values slower than the low energy 
ones. It is interesting to note that this is similar to the 
analysis of [s^ for the approach to equilibrium in plas- 
mas. 

Figure 13 provides the temporal evolution of the energy 
distribution starting from two different initial conditions. 
In both cases the system started with all particles having 
the same energy, eo = in one case, and cq = 5Tb in 
the other, and reaches its stationary distribution within 
several interactions per particle. For eo = this approach 
is uniform and faster. For eo = STb the system needs a 
slightly longer time until it arrives to the steady state, 
and the distribution behaves differently on both sides of 



eo. For e > eo single particle interactions with the bath 
dominate, and the system exhibits the bath dominated 
exponential tail. For e < gq, on the other hand, the two- 
body dissipative interactions require multiple collisions in 
order to change the shape of the distribution continuously 
from the initial delta function to the smooth steady state 
distribution. 
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FIG. 5: (Color online) Temporal evolution to the steady 
state energy distribution (dashed line) for restitution coeffi- 
cient Q = 0.5 and coupling strength / = 0.5 obtained in a 
numerical simulation of the model, starting with all particles 
having energy eo = (ascending lines after A'' and 2N inter- 
actions in the system), or eo = ^Tb (descending lines after 
iV, 2N, 3N, 4iV, and 5N interactions). 



IV. EFFECTIVE TEMPERATURES 

Various definitions of effective temperatures are used 
to characterize systems far from equilibrium. One of the 
important questions to be answered in this context is 
to what extent do effective temperatures defined by dif- 
ferent measurements on a single system yield the same 
numerical value. In Sec. Illll we first defined the granular 
temperature as the average energy per degree of free- 
dom Tg = (e), and then used the energy moments to 

(n) 

define two generalizations of Tq to the hierarchies 

and r|"'' of effective temperatures [Eq. Q]. For a non- 
exponential energy distribution these generally differ and 
moreover depend on n (see Fig. 0J|. In this section we in- 
vestigate the interrelations between three additional def- 
initions of effective temperatures - the fluctuation tem- 
perature Tp, the entropic temperature T5, and the con- 
figurational temperature Tc- 
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A. Fluctuation Temperature 

In equilibrium the fluctuation-dissipation theorem may 
be used to deduce a system's temperature from the spon- 
taneous fluctuations of any of its physical quantities. The 
theorem assures that such a fluctuation is equal to the 
temperature multiplied by the corresponding response 
function or susceptibility. As such, fluctuations and sus- 
ceptibilities of various quantities in non-equilibrium sys- 
tems may be similarly used to define effective tempera- 
tures. Here we calculate two such fluctuation tempera- 
tures and show that they generally differ one from the 
other as well as from the granular temperature Tq and 
bath temperature Tb- 

To measure fluctuation-dissipation relations, we add 
degrees of freedom {xi} to our model, whose correlations 
may be measured, and upon which a response measure- 
ment may be performed. To this end we couple the {xi} 
to an external field Fi. That is, we modify the total en- 
ergy of particle i to Ci—xiFi, and refer to {ei} as "kinetic" 
energies (or the system) and to {—XiFi{t)} as "internal" 
energies (or the probe). (This probe is added only for 
measuring fluctuation-dissipation relations as described 
in this section; all other sections of the Paper deal with 
the model defined in Sec. ITU without the "internal" de- 
grees of freedom.) 

We assume the driven dissipative dynamics defined in 
Sec. ^for the system together with non-dissipative ex- 
change of "kinetic" and "internal" energy, as described 
by Fig. ini This yields the following stochastic equations 
of motion: 



Vah 



ei{t + dt) 



e^{t) 
za [ej(t) + ej{t)] 

z {ei{t) + es] 
{e^(t) - Xi,(t)Fk{t)\ 



Probability : 
1 - Tdt 
(1 - f)Tdt (15a) 

/(I - h)rdt 

fhVdt 



Xi {t + dt) 



Value : 

Xt{t) 

z[x,{t)~ek{t)/F,{t)] 



Probability : 
1 - fhVdt (15b) 
fhVdt 



where ft, is a parameter introduced to describe the cou- 
pling strength between the system and the probe: < 
/i < 1 is the fraction of the system's interactions with 
the probe out of all its non-dissipative interactions (with 
the bath and with the probe), and k G {1, is the 
index of the particle with which particle i may inter- 
act in a "kinetic" -"internal" interaction [this is denoted 




z(e-x^F^) 
(1-z)(e.-x^F^) 



FIG. 6: The interactions between the system {e^} and the 
probe {xi}, introduced in addition to those given in Fig. Q 



differently from the index j of the second particle in a 
"kinetic" -"kinetic" interaction, since may not interact 
with Ci (thus j ^ i), while and Xi may interact {k may 
take the value i)]. 

We examine two measurements testing the relation be- 
tween steady state fluctuation and response. First, we 



consider the fluctuation (Aa; 



)^ of a sin- 



gle particle's Xi, and its response r with respect to a 
change in Fi, r = d{xi)/dFi. Second, we define the 
total system's X = ^'^'^ consider the relation 

between its fluctuation (AX^) = (X^) — (X)^, and its 
response R with respect to a change in the uniform field 
F, R = d{X) I dF. (In analogy with spin systems, Xi may 
be thought of as a single site magnetization and X as the 
total system magnetization, with (Ax^), (AX^), r and 
R the corresponding fluctuations and susceptibilities.) 

The fluctuation-dissipation theorem relates these in 
equilibrium by (Aa;^) = r • T and (AX^) = R ■ T, 
and inspires the definition in non-equilibrium systems of 
effective fluctuation temperatures Tp ^ = {Ax'^)/r and 

= (AX'^)/R for single-particle and many-particle 
measurements, respectively. In this section we calculate 
and ) for our model and demonstrate that they 
generally differ one from each other and from both Tq 
and Tb- As in we concentrate here on space- and 
time- independent fluctuation-dissipation relations; their 
temporal dependence has recently been investigated in 
j36i | and the spatial dependence in |3^ . 

By averaging Eq. lfTK|l and taking the steady state 
solution we see that 



^ N ~ F' 



(16) 



with (e) = Tb/(1 - h)/{Ai - fh) [this reduces to Eq. (gj) 
in the h ^ limit]. Therefore, 



(17) 



For the steady state averaged second moments we obtain. 



(Ax2) = ^ = 6(a,/,ft,7V)T-TG, 



2F 



(AX^) = NiAx"^) 



{XE)~{X){E) 
F 



B{aJ,h,N)-R-TG, 



(18a) 



(18b) 



where E = ^^'^ ^ B are dimensionless func- 

tions of the dimensionless model parameters and of the 
system size. We are interested in the thermodynamic 
limit {N ^ 1)]^'" which b and B reduce to the expres- 
sions given in [Sg . 

r, R, (Ax^) and (AX) all diverge as F ^ (where 
fluctuation-dissipation relations are normally measured), 
however their ratios define effective temperatures 



F 



(1) 



(Aa;2) 



6(a, /, h)Tc 



(19a) 
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^ = B{a, /, /j)rG, (19b) 

which are finite and independent of F. These fluctu- 
ation temperatures generaUy differ one from the other, 
are larger than the granular temperature Tg = (e) and 
smaller than the bath temperature Tb ■ Only in the equi- 
librium limits {a — 1 and / = 1) do all effective temper- 
atures coincide with Tb- 

is generally larger than T^^' , but in the limit stud- 
ied in [2^ of weak coupling between the system and the 
probe {h 0) the two coincide, and the expression for 
them simplifies to that given in [s^. Nonetheless, they 
differ from the granular temperature, as shown in Fig. 
[3 The difference between the two fluctuation tempera- 
tures is most prominent in the maximal dissipation limit 
(a = 0), where they reduce to the expressions given in 
|4CI|. The ratio between them for this case is shown in 
Fig. E 

The single-particle fluctuation temperature Tp"^ di- 
rectly probes the second moment of the single-particle 
energy distribution and thus gives the effective tempera- 
ture T^^^ defined in Sec. US [compare Eqs. I|16I19|I to Eq. 

(|9a|l ]. The many-particle fluctuation temperature 
on the other hand, is defined by a measurement on the 
entire system, thus reflects correlations between particle 
energies and cannot be related directly to the effective 
temperatures defined from the single-particle energy dis- 
tribution. 

In dissipative systems with strong coupling to the driv- 
ing mechanism (/ « 1) and large restitution coefficient 
(a ~ 1) the energy distribution is close to exponential, 
and correlations are weak, hence the values of Tp'' 



and Tg are similar (but not identical). As has been pre- 
dicted by kinetic theory |S|, we expect the fluctuation 
temperatures to be larger than the granular temperature 
(see Fig. [Tjl in granular gases as well, where the energy 
distribution is non-exponential. In the cases studied nu- 
merically 0, 1^ the energy distributions were only slightly 
non-exponential, resulting in small differences between 
the effective temperatures, which may explain their seem- 
ing coincidence. 
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FIG. 8: Contours of the ratio r^^'/^F ' between the many- 
particle and single-particle fluctuation temperatures vs. the 
coupling strengths / and h in the maximal dissipation limit 
{a — 0), as given in | |40|l . 



B. Entropic Temperature 
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a 

FIG. 7: Contours of the ratio Tf/Tq between the fluctuation 
temperatures and the granular temperature vs. the restitu- 
tion coefRcient a and the coupling strength / in the limit of 
weak coupling between the probe and the system (ft ^ 0), as 
given in [23. In this limit, T^^^ = T^^\ 



In analogy to equilibrium statistical mechanics, a 
further definition of an effective temperature in non- 
equilibrium systems may be constructed by differenti- 
ating the system's entropy S with respect to its aver- 
age total energy (E), yielding the entropic temperature 
Ts = {dS/d{E))-^ [illlllHlll- It is intriguing to 
inquire as to how this quantity relates to other effec- 
tive temperatures of the system and whether it has a 
fundamental thermodynamic-like significance, stemming 
from a maximization of entropy upon contact between 
systems. 

In this section we calculate Ts for the continuous en- 
ergy model defined in Sec. \^ as well as for a discrete 
energy version (described below) . We first identify simple 
scaling arguments leading to the coincidence of Tg with 
the granular temperature Tg = (e). This scaling holds 
for the continuous energy model with a > 0. We then 
demonstrate the breakdown of this scaling both in the 
singular limit of the continuous energy model at a = 0, 
which exhibits a condensation at e = 0, and by intro- 
ducing a discrete energy version of the model, where the 
structure of the particles' energetic levels contains an ad- 
ditional energy scale e. 
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1. Scaling 

For systems coupled to a heat bath of temperature Tb 
with no internal energy scale characterizing the system's 
structure (such as an interaction energy or energy spac- 
ing between possible states) the only energy scale in the 
system is Tg: When Tg is varied all energies in the sys- 
tem change linearly with its change. This is the case for 
the model defined in Sec. m Different effective tempera- 
tures, or energy characteristics of the system, may differ, 
however they all scale linearly with Tb- The ratios be- 
tween effective temperatures are dimensionless numbers 
depending on the dimensionless model parameters. 

Since the A^-particle energy distribution has dimension 
of inverse energy to the A^'th power, as long as it is non- 
singular it scales as 

(pf|J-,...,^;a,/) 
P{ei. . . . , e^) = ^ " ^, (20) 

with Lf a dimensionless function of the dimensionless en- 
ergies {ei/Ts} and the dimensionless model parameters 
[a and / in our case). 
We define the system's entropy as 

S = - J P{ei, . . . ,eAr)lnF(ei, . . . ,ejv)rfei • • ■deN-{'21) 

After changing the integration variables from {ci} to 
{ci/Tg}, and using the normalization of P(ei, . . . ,epf), 
we see that 

5 = iV In Tb + const, (22) 

with the additive constant depending only on the dimen- 
sionless parameters a and /, but not on Tb- Since the 
energy of the system {E) scales linearly with Tb [see Eq. 
Q], we may write 

S = N\n{E) + const., (23) 

and conclude that 

The functional form of the dimensionless distribution 
(p{ei/TB, . . . , gn /Tb), which depends on the model pa- 
rameters, manifests itself only in the additive term in the 
entropy, and does not affect the relation Ts — Tq . 

In order to observe richer behavior we turn to models 
with additional energy scales. We shall demonstrate this 
using a definition of Ts which is simpler to calculate, and 
which for the single energy scale case exactly coincides 
with the calculation given above. Instead of defining the 
entire system's entropy S from the entire system's en- 
ergy distribution [Eq. H21|l ]. we consider a subsystem 
comprised of a single particle and measure its entropic 



temperature. From the single-particle energy distribu- 
tion p{e), we define the single-particle entropy as 

s = — J p{e)\np{e)de, (25) 

and differentiate it with respect to the particle's aver- 
age energy: Ts = {ds/d{e))^^ . Since all particles have 
the same single-particle energy distribution, Ts is clearly 
equal for all particles in the system. The aforementioned 
scaling argument holds for the single-particle distribu- 
tion, thus for systems with a single energy scale the 
single-particle definition gives Ts — Tq as well. 

2. Condensation in the Maximal Dissipation Model 

For the maximally dissipative limit (a = 0) of our 
model, every particle undergoing an interaction with an- 
other particle in the system is left with zero energy after 
the interaction, and the system exhibits a condensation 
at e = 0. Since in the steady state for I — / of the 
particles the last interaction was such an energy drain- 
ing interaction and not an interaction with the bath, a 
fraction 1 — / of the particles have zero energy, and the 
single-particle energy distribution has the general form 

p{e) ^ {1 - me) + f r (26) 

Due to the normalization of p{e), this singular distri- 
bution yields an entropy of the form 

3 = ^(1- f)\n{l-f) + flnTB + co72st.. (27) 

The average energy per particle (e) scales linearly with 
Tb. Therefore, 

s ^ flii{e) + const., (28) 
and consequently Ts — Tq/ f, which is larger than Tq- 

3. Discrete Energy Model 

Another way to break the scaling arguments lead- 
ing to Ts — Tq is by introducing an additional en- 
ergy scale to the model. We now consider every particle 
as a quantum harmonic oscillator with possible energies 
Ci = 0, e, 2e, .... The spacing e between states constitutes 
the energy scale which invalidates the scaling arguments 
presented above. (Even in thermodynamic equilibrium, 
Ts — T IS the equilibrium temperature, while the average 
energy (e) = e/[exp(e/r) — 1] differs from the tempera- 
ture, and (e) ~ T only in the continuum limit (e) S> e.) 

The dynamics of the discrete model are as follows: The 
bath is constructed from similar quantum harmonic os- 
cillators in equilibrium at temperature Tb, thus with en- 
ergies distributed as psijie) ^ exp(— ne/Tg). In the in- 
teraction of particle i from the system with a particle 
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of energy cb from the bath their total energy ei + cb 
is conservatively redistributed between them, by ran- 
domly choosing with equal probability a new energy 
e[ S {0, e, 2e, Ci + cb}. In a dissipative interaction be- 
tween particles i and j of the system, each energy "quan- 
tum" e of the total energy ex — ei + ej has a probability 
a to remain with the interacting pair and a probability 
1 — a to be dissipated out of the system. The remaining 
energy e'j, is then randomly redistributed between the 
two particles with equal probability for every outcome 
e'i e {0,e,2e,...,e^}, and = - HH. 

We numerically solved this discrete energy model by 
Monte-Carlo simulation and obtained the average energy 
as well as the energy distribution p„ = p{ne), from which 
the entropy s = — X^^o^*" InPn was calculated. We kept 
the restitution coefficient a, the coupling strength /, and 
the energy spacing e fixed and scanned the bath temper- 
ature Tb in order to obtain the dependence of entropy on 
the average energy for given a, f and e. We compare this 
functional behavior to the corresponding equilibrium be- 
havior (i.e., conservative interactions), where the entropy 
and average energy of a quantum harmonic oscillator are 
related by 

s^^{e)) = In (^M + 1^ + M In (^^ + . (29) 

We subtracted s^'^{{e)) from the numerically obtained 
s((e)) to yield the deviation from equilibrium behavior 
of dissipative systems displayed in Fig. |51 

For (e) ^ e the discrete energy model coincides with 
the continuous energy one: For a > the entropy of the 
dissipative system in this region can be seen to merely 
be smaller by an additive constant from the entropy of 
an equilibrium system with the same energy, as expected 
from the scaling arguments for a system with a single 
energy scale [Eq. (|23|l ]. For a = 0, on the other hand, 
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FIG. 9: (Color online) The deviation of the entropy s from 
its equilibrium value s*^' vs. the average energy (e) for several 
values of the restitution coefficient a and for coupling strength 
/ = 0.5. 



the entropy deviation grows with energy, and Ts differs 
from Tq even in this high energy limit. 

We numerically differentiated s((e)) with respect to (e) 
to obtain Tg, and in Fig. 1101 compare the functional de- 
pendence of Ts on (e) to the corresponding equilibrium 
behavior. That is, for every value of (e) we normalize Ts 
by the temperature T|'' required to give this average en- 
ergy, were the system in equilibrium. Since the entropy 
of a dissipative system is smaller than that of an equi- 
librium system with the same average energy, Ts{{e)) is 
generally larger than T^'((e)). For a > 0, Tg behaves as 
in equilibrium in the two extremes of very high and very 
low energy, and exhibits non-equilibrium behavior only 
for intermediate energies ((e) ~ e). In the high energy 
limit ((e) e), the discrete energy model is equivalent to 
the continuous energy one, thus Ts behaves as in equilib- 
rium for a > 0, and reaches a value larger by a factor 1// 
from the equilibrium value for the singular limit a = 0. 

For very low energies the system behaves as a two level 
system, irrespective of whether it is in equilibrium or not. 
For the system to have such a low average energy almost 
all particles must be in the ground state (e = 0), and 
since the occupation of states rapidly decays with energy, 
only the first excited state (e — e) is relevant, while states 
of higher energy have a negligible occupation. The energy 
distribution in a two level system is characterized by a 
single number (the ratio of occupation of the two states) 
and is hence not rich enough to exhibit any features of 
a non-exponential energy distribution. For intermediate 
energies the entropic temperature exhibits significant de- 
viations from equilibrium behavior even for a > 0, as 
clearly seen in Fig. IIUI 




<e)/e 

FIG. 10: (Color online) The ratio Ts/Tj"' of the entropic 
temperature to its equilibrium value vs. the average energy 
(e) for several values of the restitution coefficient a and for 
coupling strength / = 0.5. 
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C. Configurational Temperature 

(n) 

The hierarchy of effective temperatures given by 
[Eq. H9a|) ] may be related to the hierarchy of the so-called 
hyperconfigurational temperatures. A recent extension of 
the virial theorem, states that for a system in thermo- 
dynamic equilibrium at temperature T with dynamics 
stemming from a Hamiltonian H({pi}, {qi}) (1 < « < N) 
the following relation holds 23] : 



{vn ■ B) 



(30) 



where V = {d/dpi, ...,d/dpN,d/dqi, ...,d/dqN) repre- 
sents differentiation with respect to all phase space coor- 
dinates, and B is an arbitrary vector field in phase space 
satisfying < KV"^ • < oo, < |(V ■ B)| < oo, with 
{VTi ■ B) growing slower than e^. 

If the Hamiltonian is of the form Ti. = J^iLiPi + 
V{{qi}), that is, comprised of a kinetic term depending 
only on momenta and a potential term depending only 
on coordinates, it is useful to take B = VV. This yields 
a relation between the temperature and ensemble aver- 
ages of solely the particle locations, without the need 
of measuring momenta. This may be used to define a 
configurational effective temperature in non-equilibrium 
systems. 



To = 



(|wp) 



(31) 



The hierarchy of hyperconfigurational tem- 
peratures T^"^ 1^ generalizes this by taking 
B = (0, 0, {dV/dqiy\ {dV/dqNY'), so that 



ri+1 



i=l 



(32) 



in) 

In order to interpret for our case we note that 
the model's dynamics in the non-dissipative case (a = 1) 
manifest a uniform single-particle density of states [l^ . 
That is, in the equilibrium limit our system is equivalent 
to a collection of weakly interacting harmonic oscillators, 
with ei — p1/2 + q^/2, where {qi} and {pi\ are some hid- 
den coordinates and momenta. Each qi and pi change 
periodically with time as for an isolated harmonic oscil- 
lator with energy e^. Occasionally (that is, at a frequency 
much smaller than the oscillator's frequency) this energy 
is changed due to an interaction with some other particle 
or with the bath. It is natural to extend this description 
to dissipative (a < 1) cases as well, thus we consider the 

case where V{{qi}) = Yl!i=i ^tl^- 

For a harmonic oscillator of given energy e, the tem- 
poral average of q^" over many periods of oscillation is 



e"(2n — l)!!/n!. The particle's energy changes with time 
due to interactions, thus we average over its steady state 
distribution, and obtain (q^") — (e"}(2n — l)!!/n!. Upon 
substitution in Eq. (|32|l we see that the hyperconfigu- 
rational temperatures probe ratios between succeeding 
moments of the energy distribution, and may thus be 
related to r]f^ [defined in Eq. 10]: 



T, 



Tp^"^ is undefined since averages of odd moments of q 
vanish. 

This demonstrates for one particular modeling the con- 
nection between the configurational temperature Tc^ and 
the single particle energy distribution. We suggest that 
similar relations hold in other systems, and that the hy- 
perconfigurational temperatures may be used to charac- 
terize the non-equilibrium nature of energy distributions 
in general. 



(2ri-l) 



(2n- l)(g2n-2) n(e"-i) 



(n) 



(33) 



V. DISCUSSION 

This Paper investigates several non-equilibrium phe- 
nomena observed in a minimal stochastic model for 
driven dissipative dynamics. Our model is inspired by 
granular gases, nevertheless we believe it may be rele- 
vant to a broader class of driven dissipative systems. The 
model is simple enough to admit an exact solution of the 
single-particle distribution in terms of its moments, in 
the steady state as well as during the evolution from any 
initial condition to this state. When considering particles 
with energies slightly larger than the average energy in 
the system, the high energy tail of the single-particle dis- 
tribution is seemingly overpopulated, as has been found 
in granular gases. However, we have calculated the very 
high energy tail and found that it behaves exponentially 
with a decay rate corresponding to the temperature of the 
driving bath Tq. Generally, very high energy tails man- 
ifest the bath distribution since dissipative interactions 
within the system (where typical energies are smaller 
than in the bath) hardly affect this tail. It will be in- 
teresting to investigate the tail of very high energies in 
other dissipative system, and to test whether they agree 
with the bath distributions (as has been found in lo]) 

Due to the non-exponential energy distributions of 
non-equilibrium systems, the average energy, or granular 
temperature Tg, is just one energy scale characterizing 
the system. Higher moments of the distribution define hi- 
erarchies of effective temperatures T^"'' and which 
generally vary with the order n of the moments, and co- 
incide with the actual temperature if the system is in 
equilibrium. The large n limit of these effective temper- 
atures relates to the high energy tail of the distribution. 
In our model these effective temperatures converge to the 
bath temperature in this limit, reflecting the tail's expo- 
nential behavior. Furthermore, we related these effective 
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temperatures to the hyperconfigurational temperatures 

(n) 

Tq , defined in Hamiltonian system from particle loca- 
tions. 

Fluctuation-dissipation measurements characterize as- 
pects of the second moment of the energy distribution, 
and thus show an effective temperature generally differ- 
ing from the granular temperature (or first moment). 
Here we have extended the calculations presented in 2fl| 
to general coupling strength h between the system and 
the probe, and moreover showed that single-particle and 
many-particle measurements yield different effective tem- 
peratures Tp ^ and These results, together with 
those dealing with time dependent measurements pre- 
sented in (36j , call for examining these phenomena in 
more realistic dissipative systems. 

Finally, we considered the entropic temperature Tg ob- 
tained by differentiating the entropy with respect to en- 
ergy. We showed that generally for a non-singular sys- 
tem without an internal energy scale, simple scaling ar- 
guments lead to the exact coincidence of the entropic 
temperature with the granular temperature. When this 
scaling is not valid, is generally larger than the cor- 
responding equilibrium value. It is intriguing to test this 



scaling and its breakdown in the more complex systems 
where Tg can been measured. 

We can identify an ordering of effective temperatures 
in our model: For the hierarchies defined from the en- 
ergy moments, both T^"^ and T^^^ grow with n; The 
fluctuation temperatures Tp are larger than the granular 
temperature Tg', The many-particle fluctuation temper- 
ature Tp^' is larger than the single-particle one T^^-*; All 
these effective temperatures are smaller than the bath 
temperature Tg. It will be interesting to see whether 
such ordering occurs in other driven dissipative systems 
as well. 
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